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ABSTRACT 

Observations of the magnetic fields of young solar-type stars provide a way to 
investigate the signatures of their magnetic activity and dynamos. Spectropolarime- 
try enables the study of these stellar magnetic fields and was thus employed at the 
Telescope Bernard Lyot and the Anglo-Australian Telescope to investigate two mod¬ 
erately rotating young Sun-like stars, namely HD 35296 (V119 Tau, HIP 25278) and 
HD 29615 (HIP 21632). The results indicate that both stars display rotational variation 
in chromospheric indices consistent with their spot activity, with variations indicat¬ 
ing a probable long-term cyclic period for HD 35296. Additionally, both stars have 
complex, and evolving, large-scale surface magnetic fields with a significant toroidal 
component. High levels of surface differential rotation were measured for both stars. 
For the F8V star HD 35296 a rotational shear of AH = 0.22 ;)qo 2 radd -1 was derived 
from the observed magnetic profiles. For the G3V star HD 29615 the magnetic fea¬ 
tures indicate a rotational shear of AH = 0.48 A 12 radcF 1 , while the spot features, 
with a distinctive polar spot, provide a much lower value of AH of O.OfAJg radd -1 . 
Such a significant discrepancy in shear values between spot and magnetic features 
for HD 29615 is an extreme example of the variation observed for other lower-mass 
stars. From the extensive and persistent azimuthal field observed for both targets it is 
concluded that a distributed dynamo operates in these moderately rotating Sun-like 
stars, in marked contrast to the Sun’s interface-layer dynamo. 

Key words: Line: profiles - stars: activity - stars: individual: HD 35296 - stars: 
individual: HD 29615 - stars: magnetic fields - stars: solar-type - starspots. 


1 INTRODUCTION 

Observations of the distribution of magnetic fields offer in¬ 
sights into how magnetic dynamos are governed by physical 
properties such as mass, rotation and age, with younger Sun¬ 
like stars providing proxies for understanding the history of 
the Sun’s activity and dynamo. Today we can readily ob¬ 
serve the starspots and surface magnetic fields of young Sun¬ 
like stars, to reconstruct their surface magnetic topologies 
and differential rotation. Multi-epoch observations of stellar 
magnetic regions also offers the prospect of revealing how 
magnetic activity cycles can emerge in stars like the Sun. 
Such observations can allow us to study the relative role of 
two possible dynamo theories, namely the widely accepted 
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Sun’s “shell” dynamo that operates at the interface layer 
between the radiation and convection zone and the “dis¬ 
tributed” dynamo thought to be operating in young, active 
stars where the dynamo is present throughout the convec¬ 


tion zone (e.g. Brandenburg et al. 1989 Moss et al. 1995 


Donati et al. 2003a Brown et al. 20101. While the large-scale 


toroidal magnetic field is understood to be buried deep inside 
the Sun, it has been observed on a range of rapidly rotat¬ 
ing solar-type stars through the presence of strong unipolar 


surface azimuthal magnetic holds (e.g. 

Donati et al. 

2003m 

Petit et al.j|2004a Marsden et al. 2006 

). Petit et al. 

to 

o 

o 

00 


concluded that the rotation period threshold where the sur¬ 
face toroidal field begins to dominate over the poloidal field 
is «12 d. Thus, rotation plays an important role in the gen¬ 
eration of the magnetic field and the complexity of these 
fields. 
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The study of these complex magnetic fields on Sun-like 
stars have been greatly advanced using the technique of Zee- 
man Doppler imaging (ZDI) (|Seme 1 19891 |Donati fe Semel| 


1990 Donati fe Browri||1997[ Donati et al. 2003a[ ). ZDI has 

been performed for a small sample of single, rapidly rotat¬ 
ing early G-type stars such as HD 171488 (V889 Her; Sp- 
Type: G2V) (jMarsden et al. 2006| Jeffers & Donati 2008 


Jeffers et al.|20lT ) and HD 141943 (SpType: G2) ( jMarsden 


et al. [2011a|b|. Differential rotation, given its likely impor¬ 


tance to the operation of the dynamo, is a key measurement. 
These stars have shown significant levels of differential ro¬ 
tation using both brightness and magnetic features. |Jeffers| 
& Donati (20081 observed HD 171488 and measured AH = 


0.52 ± 0.04 radd using spot features while AH = 0 .47 ± 
0.04 radd -1 was measured using magnetic features. 


Jeffers 


et al. (2011) in their more recent study measured slightly 
lower values of AH ~ 0.4 rad d -1 (Stokes /) and ~ 0.415-0.45 
rad d -1 (Stokes V). These levels of differential rotation are 
similar to an earlier study by Marsden et al. ( 2006|_of_Afi 


= 0.402 ± 0.044 radd 1 (using spot features). Marsden et 


al. (2011b) observed HD 141943 and measured AH = 0.24 
± 0.03 radd -1 using spot features with AH = 0.45 ± 0.08 
radd -1 using magnetic features on their 2010 dataset (and 
slightly lower with AH = 0.36 ± 0.09 radd -1 for their 2007 
data). Conversely, slowly rotating stars such as HD 190771 
(SpType: G2V), with a usirn of 4.3 ± 0.5 kms -1 , also ex¬ 
hibit measurable differential rota tion with AH_= 0.12 ± 0.03 
radd -1 using magnetic features ( Petit et al.|2008 ). In addi¬ 
tion to relatively high levels of differential rotation, all three 
stars exhibited complex magnetic topologies with large-scale 
poloidal and toroidal fields, each with significant higher or¬ 
der components beyond simple dipole fields. 


Barnes et al. (20051 first reported the link between 


differential rotation and effective temperature, and conse¬ 
quently convective turnover time. They were able to develop 
a correlation between rotational shear and effective temper¬ 
ature for a small sample of stars and determined that as a 
star’s effective temperature increases so does the rotational 
shear; a trend consistent with the recent theoretical calcu¬ 
lations by Kitchatinov & Olemskoy (2011). Additionally, a 
steep rise in differential rotation for late F-/early G-type 
stars appears related to the rapid shallowing of the convec¬ 


tion zone (e.g. Marsden et al. 2011bI. Kiiker, Rudiger & 


Kitchatinov (2011) demonstrate that the extreme levels of 


rotational shear observed for stars such as HD 171488 can 
only be explained with a shallow convection zone. Thus, the 
depth of the convection zone must have a major role in the 
rotational profile of a solar-type star. 

In this study, two moderately rotating, young Sun¬ 
like stars were selected based on their similar rotational 
velocities and sizes, enabling comparisons of their mag¬ 
netic fields and differential rotations. HD 35296 (SpType: 


F 8 V) is a 1. 06 ^ 0)05 M 0 (Holmberg, Nordstrom & Ander- 
son|2009 1 star with an estimated photospheric temperature 


6170K (Casagrande et al. 20111 while HD 29615 (SpType: 


G3V) is a 0.95 M s , 1.0 J? 0 star (Allende Prieto fe Lam 


bert 

1999 

). HD 35296 has 

von Eiff & Reiners 2012) w 

kms 

1 ( Waite et al. 2011b) 


of these two stars with a view to measuring the rotational 
shear and magnetic field structures. This will enable further 
comparisons with stars of similar mass but different rota¬ 


tional velocities. This is paper I of two papers focusing on 
the complex nature of the magnetic field topologies of young 
moderately rotating Sun-like stars as proxies for the evolu¬ 
tion of the young Sun. Paper II will focus on the infant Sun 
EK Draconis (Waite et al. in preparation). This paper is 
part of the BCooQ collaboration investigating the magnetic 
activity of low-mass stars (e.g. Marsden et al.|2014|. 


2 FUNDAMENTAL PARAMETERS OF THE 
STARS 

The fundamental parameters for the two targets, HD 35296 
and HD 29615, are shown in Table [T] Some of these pa¬ 
rameters such as projected rotational velocity (risim), radial 
velocity and inclination were measured as part of the imag¬ 
ing process, as explained in Sect. [5} while other parameters 
have been taken from the literature. 


2.1 HD 35296 


HD 35296 (VI119 Tau) is a F8V star ( Montes et al.|[2001[). 


It has a parallax of 69.51 ± 0.38 milliarcseconds (mas) (|van| 


Leeuwen 20 

17j) giving a distance of 14.39 ± 0.08 parsecs (pc). 

Samus et a 

. (12009) identified HD 35296 as a BY Draconis- 


type variable star. The range of age estimates of this star in 
the literature is quite large, with |Holmberg et al. (2009) esti¬ 
mating an age of 3.3 Gyr, with some authors listing this star 


as young as 20 Myr or as old as 7.5 Gyr (e.g. Chen et al. 2001 


Barry 19881. The equivalent width of the Li I 670.78 nm 
line is 99 mA, correcting for the nearby 670.744 nm Fe I 
line using the same factor developed by |Soderblom et al. 


(1993 albt. This is consistent with that measur ed by 


& Kawanomoto 


(2005) of 94.3 mA. 


Li & Hu 


Takeda 


(1998) argued 


that HD 35296 is a member of the Taurus-Auriga star form¬ 
ing region that may have reached the Zero-Age Main Se¬ 
quence (ZAMS). When placing this star on the theoretical 


isochrones of Siess, Dufour & Forestini (2000), the age of 


this star is between 20 Myr to 50 Myr. Hence it is likely 
that this star is quite youthful. 


2.2 HD 29615 


HD 29615 (HIP 21632) is a G3V star ( jTorres et al.|[2006| ). 
The HIPPARCOS space mission measured a parallax of 
18.27 ± 1.02 mas (Iva n Leeuwenl|200 7j), gi ving a distance 
of 54 . 7 / 2)9 pc. Zuckerman & Song ~ |2004 1 proposed that 
HD 29615 was a member of the Tucana/Horologium As¬ 
sociation indicating an age of ~30 Myr. Waite et al. (2011b) 


detected a magnetic field on this star, along with a varying 
emission equivalent width for the Ha line in the range from 
~370 mA to ~500 mA demonstrating the presence of a very 
active, and variable, chromosphere. 


3 OBSERVATIONS AND ANALYSIS 

HD 35296 was observed at the Telescope Bernard Lyot (TBL 
- Observatoire Pic du Midi, France) in January/February, 


1 http://bcool.ast.obs-mip.fr 
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Magnetic Fields around Moderate Rotators 3 


Table 1. The parameters used to produce the maximum-entropy image reconstructions of the two targets, including surface differential 
rotation measurements. Except otherwise indicated, parameters have been determined by this study. 


Parameter HD 35296 HD 29615 


F8V 1 G3V 2 

3.48 ± 0.01 d 3 2.34+g;g! d 4 


Spectral Type 
Equatorial Period 
Inclination Angle 

Projected Rotational Velocity, vsini 

Photospheric Temperature, T p hot 

Spot Temperature, T spo t 

Radial Velocity, v ra( j 

Stellar radius 

Stellar mass 

Age 

Convection zone depth 12 /A ( R@) 

Stokes I: f2 e q rad d 1 

Stokes 7: Af2 radd” 1 

Stokes V: f2 eq )'ari d 1 

Stokes V: Afl rad d _1 

Epoch of zero phase (MHJD) 13 


65 ± 5° 5 

cc+5 ° 

D0 -10 

15.9 ± 0.1 kms" 1 

19.5 ± 0.1 kms" 1 

6170 K 6 

5820 ± 50 K 7 

- 

3920 K 8 

38.1 ± 0.1 kms" 1 

19.33 ± 0.1 kms" 1 

1.10 R s 9 

1.0 R 0 10 

106+g;g6 M s 9 

0.95 M 0 10 

30-50 Myr 12 

30 Myr 11 

0.178 (0.201) 

0.252 (0.252) 

— 

9 co+0.06 

^• b °-0.Q2 
o n7 +0 - 10 

u ‘ u '-0.03 

1.804 ± 0.005 

n 99+0.04 

U - ZZ -0.02 

9 7/4+0.02 
-Q.Q4 

0-487°; 41 

54133.871035 (2007) 
54496.094272 (2008) 

55165.011060 (2009) 


Montes et al. 


12001 ), 


Torres et al. 


(20061, 3 Using Stokes V data, 4 Using Stokes I data, 5 The inclination angle of HD 35296 was 


based on the most recent value of P sin0 = 3.9 calculated by | Ammler-von Eiff fe Reiners| (|2012|. 6 |Casagrande et al.| |]2oTT|i, 7 Using 


the bolometric corrections of 


Bessell et al. 


(1998 \, 8 Using the relationship between photospheric and spot temperature provided by 


jBerd^ugina|j2005L 9 based on the values found in |Holmberg, Nordstrom fe Andemon|j2009 I, 10 |Allende Prieto fe Lambert] jl999| 


Zuckerman & Song (2004 1 , 12 Determined from the stellar evolution models of Siess et al.| (2000 1 , 13 Modified Heliocentric Julian Date 

(MHJD) = HJD - 24000000.0. 


2007 and again in January/February, 2008 using the high- 
resolution spectropolarimeter NARVAL. HD 29615 was ob¬ 
served at the Anglo-Australian Telescope (AAT - New South 
Wales, Australia) in November/December, 2009 using the 
high-resolution spectropolarimeter SEMPOL. Journals of 
observations are given in Tables [2] and [3] 


3.1 High Resolution Spectropolarimetric 
Observations from the TBL 

NARVAL is a bench mounted, cross-dispersed echelle spec¬ 
trograph, fibre-fed from a Cassegrain-mounted polarimeter 
unit and is similar in function to the Semel Polarimeter 
(SE MPOL) used at the AAT, which is described in Sect. 


3.2 


NARVAL has a mean pixel resolution of 1.8 kms per 
pixel with a spectral coverage from ~ 370 nm to 1048 nm 
with a resolution of -65000 spanning 40 grating orders (or¬ 
ders #22 to #61). NARVAL consists of one fixed quarter- 
wave retarder sandwiched between two rotating half-wave 
retarders and coupled to a Wollaston beamsplitter ( |Auriere1 
2003). Observations in circular polarization (Stokes V) con¬ 
sist of a sequence of four exposures. After each of the ex¬ 
posures, the half-wave Fresnel Rhomb of the polarimeter 
is rotated so as to remove any systematic effects due to: 
variations in the optical throughput, CCD inhomogeneities, 
terrestrial and stellar rotation, temporal variability and in¬ 
strumental polarization signals from the telescope and the 
polarimeter (e.g. Semel et al.||1993 Carter et al.||1996| Do 


nati et al.|2003a I. Initial reduction was completed using the 

dedicated pipeline reduction software LibreESpRIT (Echelle 
Spectra Reduction: an Interactive Tool), which is based on 
the algorithm developed by Donati et al. (1997). [Silvester 
et al. (20121 proved that NARVAL is a very stable instru¬ 


ment, with resolution and signal-to-noise ratio (SNR) being 
constant over the four years of data studies from 2006-2010. 


3.2 High Resolution Spectropolarimetric 
Observations from the AAT 


High resolution spectropolarimetric data were obtained from 
the AAT using the University College of London Echelle 


Spectrograph (UCLES) and SEMPOL (Semel et al. 1993 
Donati et al.||1997| |2003a||. The detector used was the deep 


depletion EEV2 CCD with 2048 x 4096 13.5 /im square pix¬ 
els. UCLES was used with a 31.6 gr/mm grating covering 46 
orders (orders # 84 to # 129). The central wavelength was 
522.002 nm with full wavelength coverage from 438 nm to 
681 nm. The dispersion of -0.004958 nm at order # 129 gave 
a resolution of approximately 71000. The mean resolution for 
the AAT spectra were determined tobe 1.689 km s _ 1 pixeU 1 . 
The operation of SEMPOL is similar to that explained in 
Sect. |3.ll After each exposure, the half-wave Fresnel Rhomb 
is rotated between +45° and -45° so as to remove any system¬ 
atic effects. Initial reduction was completed using ESpRIT 


developed for SEMPOL (Donati et al. 1997). SEMPOL is 


not as efficient as NARVAL, hence the SNR is lower, even 
though the AAT is a 3.9-m telescope whereas the TBL is a 
2-m telescope. In addition, the spectral range of SEMPOL 
is smaller therefore there are fewer spectral lines available 
to extract the magnetic signature from the star’s light. 


3.3 Spectropolarimetric Analysis 

Preliminary processing involved subtracting the bias level 
and using a nightly master flat combining typically 20 flat 
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Table 2. Journal of spectropolarimetric observations of HD 35296 
using the TBL. 


Table 3. Journal of spectropolarimetric observations of HD 29615 
using the AAT. 


UT Date 

UT 

Exp. Time 1 

Mean SNR for 


middle 

(sec) 

Stokes V 




LSD profiles 

2007 Jan 24 

21:08:45 

1x300 

_2 

2007 Jan 26 

21:10:23 

4x300 

19762 

2007 Jan 27 

20:18:17 

4x600 

39439 

2007 Jan 29 

19:29:15 

4x600 

38358 

2007 Feb 02 

21:17:29 

4x300 

4223 

2007 Feb 03 

20:15:21 

4x300 

29940 

2007 Feb 04 

20:06:35 

4x300 

25137 

2007 Feb 08 

20:30:51 

4x300 

25639 

2008 Jan 19 

20:01:46 

4x300 

17208 

2008 Jan 20 

21:45: 1 

4x300 

20581 

2008 Jan 21 

19:42:27 

4x300 

24563 

2008 Jan 21 

20:53:24 

4x200 

22146 

2008 Jan 22 

20:09:11 

4x300 

21913 

2008 Jan 23 

20:10:31 

4x300 

21825 

2008 Jan 24 

20:39:32 

4x300 

13296 

2008 Jan 25 

20:15:17 

4x300 

19723 

2008 Jan 26 

20:14:43 

4x300 

21148 

2008 Jan 27 

20:30:50 

4x300 

21245 

2008 Jan 28 

21:56:43 

4x300 

21245 

2008 Jan 29 

20:50:19 

4x300 

25946 

2008 Feb 02 

20:45:15 

4x300 

24375 

2008 Feb 04 

20:20:19 

4x300 

27669 

2008 Feb 05 

20:21:18 

4x300 

22909 

2008 Feb 06 

21:11: 7 

4x300 

25579 

2008 Feb 09 

20:34:17 

4x300 

24555 

2008 Feb 10 

19:13:58 

4x300 

29354 

2008 Feb 11 

20:39:11 

4x300 

25414 

2008 Feb 12 

20:45:32 

4x300 

27269 

2008 Feb 13 

20:39:30 

4x300 

24182 

2008 Feb 14 

21:00: 8 

4x300 

24386 

2008 Feb 15 

20:37: 6 

4x300 

28120 


1 Each sequence consists of four sub-exposures with each 
sub-exposure being 300 s (for example), 

2 Only one Stokes I spectrum was taken on this night. The SNR, 
as measured for order number 41, was 272. 


field exposures. Each stellar spectrum was extracted and 
wavelength calibrated against a Thorium-Argon lamp. 

After using LibreESpRIT (TBL) or ESpRIT (AAT), the 
technique of Least Squares Deconvolution (LSD) was ap¬ 
plied to the reduced spectra. LSD combines the information 
from many spectral lines to produce a single line profile, 
thereby providing an enormous multiplex gain in the SNR 


(e.g. Donati & Collier-Cameron 1997 Donati et al. 1997 


Kochukhov, Makaganiuk fe Piskunov|2010 \. A G2 line mask 

created from the Kurucz atomic database and ATLAS9 at¬ 
mospheric models ( Kurucz| |l9931 was used to compute the 
average line profile for both stars. Tables [2] and [3] lists the 
SNR for each individual Stokes V LSD profile. 

In order to correct for the minor instrumental shifts in 
wavelength space due to atmospheric temperature or pres¬ 
sure fluctuations, each spectrum was shifted to match the 
centroid of the Stokes / LSD profile of the telluric lines con¬ 
tained in the spectra, as was performed by |Donati et al.| 
(2003a I and Marsden et al. (20061. Further information on 


UT Date 

UT 

Exp. Time 1 

Mean SNR for 


middle 

(sec) 

Stokes V 

LSD profiles 

2009 Nov 25 

11:03:28 

4x750 

5567 

2009 Nov 25 

13:04:01 

4x750 

6337 

2009 Nov 25 

16:29:33 

4x750 

4908 

2009 Nov 27 

10:23:36 

4x750 

6171 

2009 Nov 27 

13:28:05 

4x750 

6978 

2009 Nov 27 

16:28:02 

4x750 

5610 

2009 Nov 28 

10:10:47 

4x750 

3821 

2009 Nov 28 

13:13:08 

4x750 

5971 

2009 Nov 28 

16:15:29 

4x750 

4280 

2009 Nov 29 

10:08:37 

400 + 3x750 

5785 

2009 Nov 29 

13:23:08 

4x750 

4604 

2009 Nov 29 

16:23:23 

4x750 

3493 

2009 Nov 30 

11:05:48 

4x750 

3283 

2009 Nov 30 

16:18:04 

4x750 

3713 

2009 Dec 01 

11:02:55 

4x750 

6353 

2009 Dec 01 

14:13:39 

4x750 

4961 

2009 Dec 02 

10:36:24 

4x750 

6158 

2009 Dec 02 

13:04:01 

4x750 

6176 

2009 Dec 02 

16:03:39 

4x750 

4651 

2009 Dec 03 

10:27:35 

4x750 

4228 

2009 Dec 03 

13:33:32 

4x750 

6312 


1 Each sequence consists of four sub-exposures with each 
sub-exposure being 750 s (for example). 


LSD can be found in Donati et al. (1997) and Kochuknov et 


al. (2010). 


4 CHROMOSPHERIC ACTIVITY 
INDICATORS 

Chromospheric activity can be determined using the Ca II 
H & K, Ca II Infrared Triplet (IRT) and Ho spectral lines. 
The two Ca II H & K absorption lines are the most widely 
used optical indicators of chromospheric activity. There have 
been a number of long-term monitoring studies of the vari¬ 
ation of Ca II H & K in solar-type stars spanning several 
decades. One study is the Mount Wilson Ca H & K project 
(e.g. Wilson 1978 Duncan et al. |1991 1. From this survey, 


differential rotation and solar-like activity cycles have been 


inferred on a number of solar-type stars (e.g. 

Baliunas et al. 

1985| Donahue & Baliunas 1992; Donahue 1993; Donahue & 

Baliunas 1994 

). 


The source functions of the Ca II H & K lines are colli- 
sionally controlled hence are very sensitive to electron den¬ 
sity and temperature. The Ca II IRT lines share the upper 
levels of the H & K transitions and are formed in the lower 
chromosphere (e.g. Montes et al.|2004l. The Ha spectral line 


is also collisionally filled in as a result of the higher tempera¬ 
tures and is formed in the middle of the chromosphere and is 


often associated with plages and prominences (e.g. Thatcher 
fe Robinson|1993 Montes et al.|2004 |. 


The Ca n H & K and Ca n Infrared Triplet (IRT) spec¬ 
tral lines were observed for HD 35296 but not for HD 29615 
as the spectral range of SEMPOL does not extend far 
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enough into the respective regions to permit monitoring of 
these diagnostic lines. 


Table 4. HD 35296: Average activity indices using the Ca II 
H&K, Ho: and Ca II IRT spectral lines. The error estimate is 
based on the range of values during each observing run. 


4.1 TBL Activity Indices for HD 35296 


Morgenthaler et al. (20121 previously reported that the spec¬ 


tra in the region of the Ca II H & K are not well normalised 
by the pipe-line reduction due to the dense distribution of 
photospheric lines. However, we did not adopt their method 
as a number of tests showed that simply removing the over¬ 
lapping section of the order provided commensurate results 
as renormalising the spectra. 

The NcanH/r-index for HD 35296 was determined using 


the method as explained in Wright et al. (2004). The result¬ 


ing Nca Hint-index was converted to match the Mount Wil¬ 
son S-values (Dun can et al. 1991) using the transformation 
shown in Eqn. [l] 


S-index = 


C1H + C2K 
C3VHK + CaRhk 


■ C s 


(1) 


where H and K is the flux determined in the line cores 
from the two triangular bandpasses with a full-width at 
half maximum (FWHM) of 0.1 nm. Two 2 nm-wide rect¬ 
angular bandpasses R hk and Vhk, centred on 400.107 and 
390.107 nm respectively, were used for the continuum flux in 
the red and blue sides of the H and K lines. The transforma¬ 


tion coefficients were determined by Marsden et al. (2014) 


by matching 94 stars that were common to both the TBL 
database and that of the Mount Wilson project. These co¬ 
efficients were: Ci = 12.873, C 2 = 2.502, C 3 = 8.877, C 4 = 
4.271 and C 5 = 1.183 x 1CT 3 . 

Two further activity indices were determined using the 
Ha spectral line and the Ca 11 IRT lines. The continuum was 
checked against that of the continuum of the synthetic nor¬ 


malised spectrum from the POLLUX database (Palacios et 
aLpOlO). In both cases,the continuum matched reasonably 
well hence was not corrected. The TBL Ha-index, based on 
Gizis, Reid & Hawley (2002), was determined using Eqn. [ 2 ] 


Index 

2007 

2008 

S-index 1 

0.323 ± 0.008 

0.322 ± 0.016 

N He, 

0.310 ± 0.004 

0.314 ± 0.004 

N CaIRT 

0.887 ± 0.011 

0.891 ± 0.015 


1 The resulting N ( y, J n//^-index was converted to match the 
Mount Wilson S-values (D uncan et al.j 1 991| using the transfor¬ 
mation shown in Eqn. [T] 


(right panels) for HD 35296. The epoch for the 2007 data was 
set to the middle of the observation run whereas the epoch 
of the 2008 data was set to 104 stellar rotations later on and 
was chosen to be close to the middle of the 2008 observation 
run, as explained in Sect. |5.2| The index was measured for 
each spectrum, then combined in sets of four, coinciding 
with a cycle of 4 sub-exposures as explained in Sect. |3.l] The 
average of these four indices was determined with the error 
bar being the minimum and maximum values for that set. 
The data from 2007 showed limited variation, either due to 
the reduced coverage or a more homogeneous chromosphere. 
In contrast, in 2008, HD 35296 exhibited greater modulation 
of the S-index, CaIRT-index and Ha-index while the average 
remained very similar with an S-index for 2007 of 0.323 ± 
0.008 and for 2008 of 0.322 ± 0.016. This is consistent, within 
error limits, with the average S-index of 0.308 ± 0.07 from 


the Mount Wilson survey (Duncan et al. 

1991). Using the 

transformation values from Rutten 

(1984 

1 , Log R' hk was 


determined to be -4.45 ± 0.03. This was an average value 
of the S-index over the observing runs with the error based 
on the minimum and maximum values. Table [4] show the 
average S-index, N Ha and N caiRT indices for 2007 and 2008. 
Each value is the average over the observing run with the 
error showing the variation during this time, most likely as 
a result of modulation due to the rotation of the star. 


NHa-index = 


Fhc 


( 2 ) 


ITfa "t Rlla 

where F Ha is the flux determined in the line cores from 
the triangular bandpass with a FWHM of 0.2 nm. Two 
0.22 nm-wide rectangular bandpasses Vhc, and Rrr a , cen¬ 
tred on 655.885 and 656.730 nm respectively, were used for 
the continuum flux in the red and blue sides of the Ha line. 
The TBL CaIRT-index, following Petit et al. (20131, was 
determined using using Eqn. [3]. 


NcaiRT-index = — —— (3) 

VlRT + ttlRT 

where Y, Firt is the total flux measured in the line cores of 
the three spectral lines, 849.8023, 854.2091 and 866.2141 nm 
using triangular bandpasses with a FWHM of 0.1 nm. Two 
0.5 nm-wide rectangular bandpasses Virt and R/kt, cen¬ 
tred on 847.58 and 870.49 nm respectively, were used for the 
continuum flux in the red and blue sides of the Ca 11 IRT 
spectral lines. 

Fig .[TJshows the variation in the Ca 11 H & K lines, Ca 11 
IRT and Ha spectral lines during 2007 (left panels) and 2008 


4.2 AAT Ha activity index for HD 29615 


The Ha spectral line was used as a chromospheric indica¬ 
tor for HD 29615. A similar approach, as explained in |4.1| 
was applied to the AAT data. The overlapping sections of 
the AAT spectra were removed and the continuum checked 
against a synthetic spectrum obtained from the POLLUX 
database, In all cases, the continuum near the Ha spec¬ 
tral line closely matched the synthetic spectra. Eqn. [2] was 
then used to determine the Ha index for HD 29615. Fig. [2] 
shows how this Ha-index varied over the rotation period of 
HD 29615. The error bars indicate the range of the mea¬ 
surement during a cycle of 4 sub-exposures, as explained in 
Sect. |3.1~| This supports the variable nature of the mid-level 


chromosphere, first observed by Waite et al. (2011b). 


5 IMAGE RECONSTRUCTION 

Surface images of these moderately rotating stars were gen¬ 
erated through the inversion of a time series of LSD pro¬ 
files of Stokes I (brightness images) or Stokes V (magnetic 
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2007 


2008 



Figure 1. The variation in the Ca II H & K lines (upper), Ha (middle) and Ca II IRT (lower) spectral lines during 2007 (left) and 2008 
(right) for HD 35296. The index was measured for each spectrum, then combined in sets of four, coinciding with a cycle of 4 sub-exposures 
as explained in Sect. The average of these four indices was determined with the error bar being the minimum and maximum values 
for that set. The rotational period used was 3.48 d. The observing run covered ~4.3 rotations in 2007 and ~7.7 rotations in 2008. The 
symbols represent the respective rotations: . only = first, o = second, □ = third, A = fourth, v = fifth, <3 = sixth, [> = seventh and ☆ 
= eighth rotation. 


images). The Zeeman Doppler imaging code used was that 
of Brown et al. (19911 and Donati & Brown (19971. Be¬ 
cause this inversion process is an ill-posed problem where 
an infinite number of solutions is possible when fitting to 
the noise level, this code implements the Skilling fe Bryan] 
maximum-entropy optimization that produces an im¬ 


age with the minimum amount of information required to fit 
the data to the noise level. 


5.1 Brightness image of HD 29615 

Brightness images were not reconstructed for HD 35296 
as the risim of 15.9 kms" 1 is below the usual limit for 
Doppler imaging. In contrast, the observed intensity pro¬ 
files for HD 29615 contained sufficient information for both 
mapping and a measurement of differential rotation using 
spot features. 

A two-temperature model, one being the temperature 
of the quiet photosphere while the second is that of the cool 
spots, has been used to reconstruct the brightness image for 
HD 29615. Synthetic Gaussian profiles were used to repre- 
the spot and photosphere. |Unruh| 
) showed that the use of these syn¬ 
thetic profiles provides brightness maps that are very similar 


sent the profiles of both 


& Collier Cameron (1995 


to those maps created by using profiles from slowly rotating 
stars of commensurate temperature. Hence, many authors 
now adopt this approach (e.g. |Petit et alT] 


|Ma 


den et al. 2005| 2006 2011a|b Waite et al.| 2011a). Using 
the relationship between photospheric and spot tempera¬ 
ture provided by |Berdyugin a ([2005]), the spot temperature 
of HD 29615 was estimated to be 3920 K. 

The imaging code was used to establish the values of a 
number of basic parameters, including the star’s projected 
rotational velocity, vsini, and radial velocity, v ra d- This was 
achieved by systematically varying each param eter in or der 
to minimize the reduced-y 2 (Xr) value (e.g. Marsden et al. 
2005 Jeffers & Donati 2008). These key stellar parameters 


were determined with initially v ra£ j and then Dsini This se¬ 
quence was repeated each time additional parameters, such 
as inclination or differential rotation parameters, were mod¬ 
ified as a result of the imaging process. 

The inclination angle for HD 29615 was estimated us¬ 
ing the bolometric corrections of |Bessell, Castelli fe Plez~| 
(1998) and was confirmed using the imaging code. This value 
was determined to be 65|o°. The minimum angle of ~ 50° 


assumed a stellar radius of 1.2 R 0 (Messina et al. 2010). 


A lin ear limb-darkening coefficient of 0.62 was used (Sing 
2010). The full set of parameters that gave the minimum Xr 
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Figure 2. The top panel represents the variation in the Ha in¬ 
dex as a function of longitude for HD 29615 (see Seel.. [4~~2] l. The 
observing run covered ~3.5 rotations. The symbols represent the 
respective rotations: . = first, o = second, □ = third and A = 
fourth rotation. The bottom three panels represent the variation 
in magnetic field strength for each configuration as a function of 
longitude (see Sect. |6. 1 [ |. The shaded regions show the variation 
in the strength at each longitude as determined by varying usini, 
inclination angle, Q eq and Af2 within the limits of their respec¬ 
tive accuracy. The Epoch was set to MHJD = 55165.011060 with 
a rotational period of 2.34 d. 


value of 0.45 are shown in Table [T] and were adopted when 
producing the final Doppler imaging map as shown in the 


JJopp 

top left panel of Fig. Isl A \r value less than one is pos¬ 
sible as the SNR calculated for the LSD Stokes I profiles 
are underestimated (e.g. Petit et al. 2004b; Marsden et al. 


201 la I. Nevertheless, this has no impact on the final bright¬ 


ness maps. The maximum-entropy fits to the Stokes I LSD 
profiles for HD 29615 with the measured surface differen¬ 
tial rotation incorporated into the analysis is shown in Fig. 
[4] Deviations of the model profile from the observed profile 
are shown in the dynamic spectrum displayed in Fig. [5] To 
produce this dynamic spectrum, each observed LSD profile 
was normalised, then subtracted from the associated nor¬ 
malised model profile with the residual profile being plotted 
as the dynamic spectrum. The scale of this dynamic spec¬ 
trum is ± 4 x 1(T 3 or ss 0.4 per cent of the intensity of the 
line profile. The small residuals indicate that the majority 
of the large-scale features can be accurately modelled. How¬ 
ever, there is a consistent mis-match between the modelled 
and original data as evidenced by a brighter band near the 
red wing of the dynamic spectrum. This indicates that the 
modelled data has not incorporated information from the 
wings of the original line profiles. This is advantageous in 
that modelling of that part of the line profile may have led 


to spurious banding in the resulting image (Unruh & Collier 
Cameron||1995 1. 
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5.2 Magnetic images 

The magnetic topology was determined for both HD 35296 
and HD 29615 using ZDI. ZDI requires high SNR data as 
the polarization signature is typically less than 0.1 per cent 
of the total light intensity (D onati et al.( l997). The Stokes 
V data was used to reconstruct radial, azimuthal and merid¬ 
ional fields coinciding with the three vector fields in spherical 


geometry. The modelling strategy of Donati & Brown (19971 


was used to construct the magnetic field topology on both 
stars. ZDI only measures the large-scale magnetic field as 
the small-scale magnetic fields cannot be recovered as the 
positive and negative magnetic fields within the resolution 
element are likely to cancel each other out. The mapping 
procedure involved uses the spherical harmonic expansions 
of the surface magnetic field, as implemented by |Donati et al.| 
12006). The maximum spherical harmonic expansion imax 


= 11 was selected for HD 35296 while imax = 14 was used 
for HD 29615. These were the minimum values where any 
further increase did not produce any difference in the mag¬ 
nitude and topology of the magnetic field recovered. 

The magnetic maps for HD 29615 are shown in Fig. [3] 
while the associated fits between the modelled data and the 
actual LSD profiles are shown in Fig. [g] The Xr value of the 
magnetic models for HD 29615 was 1.25; showing that the 
fit accuracy was very close to reaching the noise level of the 
data. 

The magnetic maps for HD 35296 for the 2007 and 2008 
datasets are shown in Fig.[7]with the associated fits between 
the modelled data and the actual profiles shown in Fig. [8] 
The LSD profile taken on February 02, 2007 of HD 35296 
was excluded from the mapping process due to the rela¬ 
tively poor SNR. Including these data made no difference 
to the resulting magnetic maps and field configuration but 
did make a marginal degradation of the total magnetic field 
strength (by 0.1 G). In each case, the minimum \r value of 
the magnetic models was set to 1.0. 

The full set of parameters used to produce these maps, 


including differential rotation (see Sect 5.31 are listed in Ta¬ 


ble [T] Due to the length of time between the observing runs 
for HD 35296, the phases for the 2007 run were calculated 
using an epoch of the Modified Heliocentric Julian Date 
(where MHJD = HJD - 2400000.0) of 54133.871035 whereas 
the 2008 run used an epoch of MHJD = 54496.094272. This 
2008 epoch was selected as it was close to the mid-point of 
the observing run, but was an integer number of rotations 
since the mid-point of the 2007 run based on the 3.48 day 
period. 


5.3 Differential rotation 


There have been a number of different methods by which 
differential rotation has been measured on stars. IReiners ~fcl 
Schmitt (200 3a|b ) have used a Fourier transform method 
to derive the parameters for rapidly rotating inactive F- 
type / early G-type stars. They produce a deconvolved line 
profile from the stellar spectrum and determine the ratio 
of the second and first zero of the resulting Fourier trans¬ 
form. This ratio is a measure of the magnitude of the dif¬ 
ferential rotation on the stellar surface. However, for active 
young stars with significant asymmetry within the line pro¬ 
files, this technique is not as effective. Differential rotation 
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Figure 3. The maximum-entropy brightness and magnetic image reconstructions for HD 29615. These maps are polar projections 
extending down to -30°. The bold lines denote the equator and the dashed lines are +30° and +60° latitude parallels. The radial ticks 
indicate the phases at which this star was observed and the scale of the magnetic images is in Gauss. Differential rotation, as measured 
using the Stokes I and V profiles respectively, has been incorporated into the reconstruction of the maps. The spot map has a spot 
coverage of 2.6 per cent while the global magnetic field strength is 81.6 G. The Epoch was set to MHJD = 55165.011060 with a rotational 
period of 2.34 d. 


of active, rapidly-rotating stars can be estimated by track¬ 
ing spot features at different latitudes on the surface of the 


star from multiple Doppler images (Collier Cameron, Donati 


& Semcl 2002). Alternatively, cross-correlation techniques 


can be used on two independent Doppler images (Donati & 
Collier-Cameron] l9971. Finally, a solar-like differential rota¬ 
tion law, as defined in Eqn. [fj can be applied to measure the 
rotation of a number of solar-type stars (e.g. Petit, Donati 


& Collier Cameron|2002; 

Donati, Collier Cameron & Petit 

2003b Petit et al. 2004b 

Barnes et al. 2005| Marsden et al. 


2005 2006). 


0.(8) = O eq - A Osin 8 


(4) 


where 12 ( 0 ) is the rotation rate at latitude 8 in radd -1 , 12 eq 
is equatorial rotation rate, and Afl is the rotational shear 
between the equator and the pole. This technique uses fixed 
information content of the spot coverage (Stokes I) or mag¬ 
netic field strength (Stokes V) for the maps while systemati¬ 
cally adjusting the differential rotation parameters, £ 2 eq and 
AO. The best fit to the data, using \ 2 minimization tech¬ 
niques, is subsequently determined by fitting a paraboloid to 
the reduced-x 2 landscape thereby measuring the differential 


rotation parameters. It is this x ' 2 minimization technique 
that has been adopted for this work as it is best suited to 
the longer time-base data sets presented in this paper. 


Using the magnetic signatures, HD 35296 has an equa¬ 
torial rotational velocity, O eq , of 1.804 ± 0.005 radd -1 with 
a rotational shear, AO, of 0.22*o;o| radd -1 . Using the bright¬ 
ness features, HD 29615 has an equatorial rotational veloc¬ 
ity, o eq , of 2.6818:85 radd 1 with a rotational shear, AO, 
of 0.07lo;o3 radd -1 . Using the magnetic features, HD 29615 
has an equatorial rotational velocity, O eq , of 2.74lo;o 2 radd -1 
with a rotational shear, All, of O. 48 I 01 1 radd -1 (see Table 

0 


The errorsfor the differential rotation measurements 
were calculated by individually varying stellar parameters 
such as spot occupancy (for Stokes /: ± 10 per cent) or global 
magnetic field (for Stokes V: ± 10 per cent), inclination an¬ 
gle (± 5° for HD 35296 and +5 and -10 ° for HD 29615) and 


risini (± 0.1 kms -1 ) and determining the minimum f 2 e 9 -A 12 

g air f rom each of the reduced \ 2 landscapes generated. Figs. 
and[ 


10 


show the reduced \ 2 landscape for the optimum set 
of parameters, as listed in Table [T] Superimposed on these 
are the individual O eq -AO pairs derived from this variation 
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Figure 4. The maximum-entropy fits to the Stokes / LSD profiles for HD 29615 with the measured surface differential rotation 
incorporated into the analysis. The red lines represent the modelled lines produced by the Doppler imaging process whereas the black lines 
represent the actual observed LSD profiles. Each successive profile has been shifted down by 0.020 for graphical purposes. The rotational 
cycle at which the observation took place are indicated to the right of each profile. The Epoch was set to MHJD = 55165.011060 with a 
rotational period of 2.34 d. 


of stellar parameters, with the error bars being 1-cr errors 
in the paraboloid fit. The final ellipse was generated to en¬ 
compass all the differential rotation values. The true error 
bars may be slightly larger than this due to intrinsic spot 


evolution during the data collection (e.g. Morgenthaler et 
|aL||2012). The differential rotation parameters determined 
for HD 35296 for the 2008 observing run were used in the 
mapping process for the 2007 data as no differential rotation 
was able to be determined for this dataset due to the limited 
number of exposures. 


6 DISCUSSION 

HD 35296 and HD 29615 are two moderately rapidly rotat¬ 
ing young Sun-like stars that display rotational variation in 
their chromospheric activity. In this paper, we have deter¬ 
mined the surface topology of both stars revealing complex 
magnetic fields and high levels of surface differential rota¬ 
tion. 


6.1 Chromospheric activity 

Chromospheric activity was determined using the Ca II H 
& K, Ca II Infrared Triplet (IRT) and Ha spectral lines for 
HD 35296 while only the Ha line was used for HD 29615 as 
the spectral range of SEMPOL did not extend far enough 
into the respective regions to permit monitoring of the Ca II 
diagnostic lines. Table [4] shows the average activity indices 
in 2007 and 2008 for HD 35296. The average S-index was 
0.322±0.016 while the average Ha index, Nrf a , 0.312 ± 0.006 
and the average Ca II IRT index, Nc 0 irt , was 0.889lo;oi3. 
Fig.g demonstrates the variable nature of the activity of 
HD 35296. In 2007, the chromosphere was more homoge¬ 
neous with respect to phase. The large error bar at phase ~ 
0.15 is the result of a poor signal-to-noise observation. The 
amplitude of the modulation in the Ca II H & I< and Ha 
indices is less than that observed in 2008. This would indi¬ 
cate that HD 35296 is more active during the 2008 epoch. 
In addition, the error bars on each of the data points do 
not overlap indicating that the chromosphere not only var¬ 
ied as a function of rotation, but also from one rotation to 
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Figure 5. Dynamic spectrum of the residual plot for the Stokes 
I LSD profiles for HD 29615. Each LSD profile was normalised 
and subtracted from the associated normalised model profile. 



the next. This is clearly observed in all three chromospheric 
diagnostic lines. The conclusion is that the chromosphere of 
HD 35296 is highly variable over a one month interval of 
approximately eight stellar rotations. 

S-index measurements are also available from IBaliunasI 


et al. (19951, where data was obtained between 1966 to 1991. 


The average S-index from the Baliunas sample is 0.33 ± 0.02 
and the Log R’hk is -4.36 ± 0.04. Mount Wilson data is not 
available from 1991 onwards but fortunately Lowell observa¬ 
tory started a long-term Ca II H & K monitoring programme 


where HD 35296 was observed from 1994 to 2008 (Hall et 


al. 2009). All available S-index values scaled to the Mount 


Wilson system are shown in Fig. |11| HD 35296 exhibits a 
gradual decrease in its S-index from 1966 to 1980 and shows 
a relatively flat trend from 1980 onwards. Interestingly dur¬ 
ing the low activity period the S-index does exhibit small 
variations. The long-term variations also indicate a proba¬ 
ble long cyclic period which was also determined by |Baliunas| 
|et al.|fl995| ). No error estimates are available for the Lowell 
catalogue. 

The average Ha index for HD 29615 is 0.373 ± 0.008. 
As indicated in Sect. |4j the Ha spectral line is often associ¬ 


ated with plages and prominences (e.g. Thatcher & Robin- 
1993 Montes et al. |2004 1. There was evidence of the 


existence of plage-like behaviour on HD 29615 when con¬ 
sidering the Ha emission as shown in Fig. [2] (top panel). If 
longitude increases with phase with 0° longitude correspond¬ 
ing to 0.0 phase, then this Ha emission reaches a maximum 
at longitude ~220°. There appears to be enhanced magnetic 
features in the radial magnetic field trailing slightly, at ap¬ 


proximately 180°, as shown in Fig. [2] The three lower pan¬ 
els in Fig. [2] show the fractional magnetic field strength as 
a function of longitude around the star for the radial, az¬ 
imuthal and meridional magnetic field components. These 
were calculated using Eqn. [5] 


F(tp) = B(tp)dtp (5) 

where F(<p) is the fractional magnetic field strength at lon¬ 
gitude 1 p, B(<p) is the average magnetic field strength at lon¬ 
gitude ip, and dp is the width of each latitude ring. These 
enhanced features can also be seen in the resulting map in 
Fig. [ 3 ] (top right panel). We can speculate that the enhanced 
positive radial magnetic field is contributing to the activity 
in the mid-level chromosphere, perhaps with the production 
of plages or even prominences at mid-latitude, very similar 
in location to those features observed on the Sun. 


6.2 Brightness and Magnetic Maps 


The results presented in this paper show that HD 35296 
and HD 29615 are young, Sun-likc stars whose moderately 
rapid rotation has led to very complex large-scale surface 
magnetic field topology. As discussed in Sect. [5} brightness 
images were not produced for HD 35296 due to its relatively 
low usim with the LSD intensity profiles not being suffi¬ 
ciently deformed by any surface features present to recover 
reliable spot information using DI. The brightness map for 
HD 29615, in Fig. [ 3 ] (top left panel), shows a dominant po¬ 
lar spot. This polar spot is slightly off-centred with features 
extending to ~60° latitude for phases between, </>, 0.0 to 0.5, 
but not evident for (f> = 0.5 to 1.0. While the polar spot is 
very prominent, there appears limited spot coverage at mid- 
to lower-latitudes with a small number of features observed 
at (j) ~ 0, 0.25-0.3 and 0.55-0.65. The features at (ji ~ 0 and 
~ 0.25 may not be reliable as the dynamic spectrum shown 
in Fig. [5] displays a dark diagonal band moving from the 
blue wing to the red wing. This indicates that the modelled 
data may have incorporated spots at this phase that were 
not present in the original data hence could not accurately 
constrain the latitude of these two features. lUnruh fe Collier] 


Cameron (1997) note that spectroscopy extracts the higher 


latitude features but cannot discriminate lower latitude fea¬ 
tures at latitudes < 30°. However, the paucity of spots at 
mid- to lower-latitudes is dissimilar to other more rapidly 
rotating solar-type stars like He 699 (V532 Per, SpType: G2- 
3V) ( Jeffers, Barnes fe Collier Cameron||2002 1, HD 106506 
(SpType: G1V) (jWaite et al. 2011a I and HD 141943 (Sp¬ 
Type: G2) (Marsdc n et al.|2011a| ) that appear to have spots 
covering all latitudes and phases. 

The magnetic maps for HD 29615, in Fig. [3j show az¬ 
imuthal, meridional and radial magnetic field structures. 
The azimuthal magnetic field is strongly positive although 
not a complete ring of field, unlike that observed on other 
solar-type stars such as HD 106506 ( Waite et al.|2011a I and 
HD 141943 (Marsden et al. 2011al. The radial field is also 


strongly positive at higher latitudes, from approximately 
+30° to the pole. The maps for HD 35296 from 2007 and 
2008 are shown in Fig. [7] Like HD 29615, HD 35296 dis¬ 
plays a predominantly positive azimuthal magnetic field in 
both 2007 and 2008 while the radial and meridional fields 
show mixed polarity across the surface. 
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Figure 6. The maximum-entropy fits to the Stokes V LSD profiles for HD 29615 with the measured surface differential rotation 
incorporated into the analysis. The black lines represent the observed Zeeman signatures, while the red line represent the modelled lines. 
Each successive profile has been shifted down by 0.002 for graphical purposes. The minimum Xr value was set to 1.25. The rotational 
phases at which the observations took place are indicated to the right of each profile. The error bars to the left of each profile are ± 0.5 
(j. The Epoch was set to MHJD = 55165.011060 with a rotational period of 2.34 d. 


When using Stokes V without Stokes /, Q, U there is 
crosstalk between the three magnetic field configurations, 
although the azimuthal field is the most robust with lim¬ 
ited crosstalk. At high latitudes, ZDI can recover all three 
field components but as lower latitudes, in particular for 
stars with inclination angles greater than 50°, the crosstalk 


is from the meridional to radial field component (Donati & 
Brown|1997 1. This crosstalk appears stronger for HD 35296 


when compared with HD 29615. This could be due to either 
a weaker magnetic field or fewer phases for each rotation for 
HD 35296 thereby affecting the code’s ability to distinguish 
between the radial and meridional field components, partic¬ 
ularly at lower latitudes. Rosen & Kochukhov (2012) have 
observed when using all four Stokes IQUV parameters, the 
crosstalk between the components is significantly reduced. 


However, the linear polarization state (Stokes Q and U pa¬ 
rameters) are significantly weaker than the Stokes V and are 
extremely difficult to observe on active cools stars such as 
HD 29615 and HD 35296 (e.g. Wade et al.|2000 Kochukhov 
et al.|2004 1. 


6.3 Latitude dependence of the magnetic fields 

The fractional magnetic field strength as a function of lati¬ 
tude for each orientation for HD 35296 and HD 29615 was 
determined using Eqn. [6] 


F(6) 


B(8)cos(8)d9 


( 6 ) 
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Figure 7. Magnetic imaging of the large-scale field for HD 35296 from both the 2007 (top series) and 2008 (bottom series) observing 
runs. These maps are polar projections, as explained in Fig. [3] The 2007 global magnetic field was 13.4 G whereas the 2008 field was 
17.9 G. The Epoch was set to MHJD = 54133.871035 for the 2007 data and MHJD = 54496.094272 for the 2008 data, with a rotational 
period of 3.48 d. 


TBL: 2007 TBL: Jan 2008 TBL: Feb 2008 
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Figure 8. The maximum-entropy fits to the LSD profiles for HD 35296 during 2007 (left panel) and 2008 (right panel), with differential 
rotation incorporated into the analysis. The observed Zeeman signatures are shown in black while the fit to the data is shown as red lines. 
The rotational cycle and ± 0.5 a error bars of each observation are shown next to each profile. Each profile has been shifted down by 
0.0012 for clarity. The minimum \r value used was 1.0. The LSD profile obtained on February 02, 2007 was excluded from the mapping 
process due to its relatively low signal-to-noise. The rotational period used was 3.48 d. 
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Figure 9. Differential rotation using Stokes V: HD 35296 in 
2008. The contour plot is a 1-cr projection using the optimum 
set of parameters, as listed in Table [T] The darker regions corre¬ 
spond to lower reduced-x 2 values. Superimposed on this grid is 
the error ellipse (thick line) that encompasses a range of n eq - 
An pairs determined by varying such parameters as vs'mi, mag¬ 
netic field strength and inclination. Each datapoint is the result 
of paraboloid fit of the reduced-x 2 landscape for each individual 
differential rotation value. The error bar on each datapoint are 
1-cr errors in the paraboloid fit. 
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Figure 10. Differential rotation using Stokes V: HD 29615. The 
contour plot is a 1-cr projection using the optimum set of param¬ 
eters, as listed in Table [I] Superimposed on this grid is the error 
ellipse (thick line) encompassing the respective n eq - Af2 pairs, 
as explained in Fig. [9] The error bar on each datapoint are 1-cr 
errors in the paraboloid fit. 


where F(0) is the fractional magnetic field strength at lati¬ 
tude 6 , B (0) is the average magnetic field strength at lati¬ 
tude 0 , and d6 is the width of each latitude ring. The results 
are shown in Fig. |12[ In order to determine the robustness 
of each measurement, the vsi m, v ra d,inclination angle, Q eq 
and were varied, placing limits on the fractional mag¬ 
netic field at each latitude for each of the orientations. These 
measured variations are shown by the shaded regions while 
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Figure 11. The long-term Ca II H & K S-index data for 
HD 35296. The Mount Wilson data taken from 1966 to 1991, 
shown as • j jBaliunas et al.]|1995[) while th e Lowell data taken 
from 1994 to 2008, shown as t ( [Hall et ah]|2009[ >. Additionally, 
the two values shown as m were determined by this work. 
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the values for the optimum parameter set are shown as a 
solid line for the 2007 data and a dotted line for the 2008 
data. HD 35296 displayed both positive and negative ra¬ 
dial fields in 2007 but predominantly positive held in 2008 
whereas the azimuthal held remains a positive held for both 
epochs. Conversely, the radial held, at mid-latitudes, ap¬ 
pears to have undergone a reversal of polarity from 2007 to 
2008. This may be real or an artefact due to fewer prohles 
being used in the reconstruction. 

To test this, we reduced the number of prohles in the 
2008 dataset, to match the same number of prohles and 
phase as that obtained in 2007. This had an impact on 
the strength of the global magnetic Held (17.9 G reduces 
to 11.9 G) and the latitudinal distribution of the respec¬ 
tive Helds. The azimuthal held remained strongly positive, 
as shown in Fig. 12 (dot-dash line with cyan colour). How¬ 
ever, the radial held no longer shows a strong positive peak 
at ~ 70° latitude when compared with using the full number 
of prohles, although it does remain positive at these mid¬ 
latitudes. As a further test, maps were produced from the 
hrst two rotations, middle four rotations and final two rota¬ 
tions. This demonstrated the robust nature of the azimuthal 
magnetic held, although fewer prohles meant that the imag¬ 
ing code could not adequately recover the latitude of each 
feature. Nevertheless, the mid-latitude reversal in the ra¬ 
dial held from 2007 to 2008 may well be real. The evidence 
is tantalising but more data would be required before any 
dehnitive conclusion can be drawn regarding any changes in 
the mid-latitude held polarity of HD 35296. 

Comparing HD 35296 with HD 29615, it appears that 
the dominant azimuthal Helds are restricted to the equato¬ 
rial regions for HD 35296 whereas the similar Helds domi¬ 
nate at higher latitudes for HD 29615. Both HD 35296 and 
HD 29615 have similar latitudinal distributions in the merid¬ 
ional and radial Helds. Thus both of these young, Sun-like 
stars with similar masses and t;sini have similar magnetic 
held configurations. 


6.4 Magnetic field configurations 

Large-scale magnetic Helds on the Sun are considered to arise 
from the interplay between the poloidal and toroidal com¬ 
ponents in what is commonly referred to as the a - ff dy- 
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Figure 12. Variation of the respective magnetic field orientations for HD 35296 (left series of panels) and HD 29615 (right series of 
panels). The shaded regions show the variation in the strength at each latitude point as determined by varying usini, inclination angle, 
fl eq and AD within the limits of their respective accuracy (see Sect. |5.3| l. Considering HD 35296, the solid line with the dark blue shading 
is for 2008 data while the dotted line with the green shading is the 2007 data. The dot-dash line with cyan shading represents the 2008 
data with reduced number of profiles approximately matching the number of profiles and phase of the 2007 data. The period used for 
HD 35296 was 3.48 d while the period used for HD 29615 was 2.34 d. 


namo (Babcock 1961). In the Sun, the toroidal component 
is understood to be confined to the interface between the 
radiative core and the convection zone (e.g. |MacGregor fe| 


Charbonneau 1997 Charbonneau & MacGregor 19971. How¬ 


ever, in solar-type stars, the toroidal component may man¬ 
ifest itself in the form of strong azimuthal magnetic fields 


on, or near, the surface of the star (e.g. Donati et al.|2003a 


Petit et al. 2004b I. Additionally, Petit et al. (20081 infer 


that a rotation period lower than ~12 d is necessary for 
the toroidal magnetic energy to dominate over the poloidal 
component. Both of these moderately rapidly rotating stars 
appear to have a large azimuthal component. These results 
are consistent with observations of more rapidly rotating 
stars exhibiting similarly strong azimuthal magnetic fields 
such as AB Doradus (SpType: K0V), LQ Hydrae (SpType: 


K0V) and HR 1099 (SpType: K2:Vnk) (Donati et al. 2003a), 

HD 171488 (SpType: G2V) ( 

Marsden et al. 2006; Jeffers 

& Donati||2008j Jeffers et al. 

2011), HD 141' 

M3 (SpType: 
pType: G1V) 
trong surface 
e been inter- 

G2) (iMarsden et al. 

2011 b) an 

d HD 106506 (S 
tions of these s 
active stars hav 

(Waite et al. 2011a 
azimuthal magnetic 

. Observa 
fields on 


preted by Donati et al. (2003a I as a result of the underlying 
dynamo processes being capable of generating fields directly 
in the subsurface region (e.g. |Dikpati~et al. 2002) and could 
be distributed throughout the whole convection zone (e.g. 
Lanza, Rodono &: Rosner]l9981. 


The magnetic field on HD 35296 was predominantly 
poloidal in 2007, with ~82li per cent of the total magnetic 
energy. By 2008, the field had reorganised itself to a ~50:50 
per cent poloidal-toroidal configuration. This is shown in 
Table [5] and in Fig. [7] Again there are two plausible rea¬ 
sons for this. One is that the magnetic topologies had sig¬ 
nificantly changed. Alternatively, simply using fewer pro¬ 
files may be responsible for the observed changes. Reducing 
the number of profiles, as explained in Sect. |6.3| had only 
a minor effect on the balance of the poloidal-toroidal field 
configuration in 2008. This is shown in Table [5j where the 
full dataset produced a poloidal-toroidal ratio of 50:50 per 
cent while the reduced dataset produced a ratio of 48 to 
52 per cent. Higher cadence data would be required to sup¬ 
port the conclusion that the global magnetic field topology 
had significantly evolved during the course of one year on 
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Table 5. Magnetic quantities derived from the set of magnetic maps for HD 35296 and HD 29615. Both the poloidal and toroidal 
components are listed along with the relative weightings in each geometry. The errors are variation bars, as described in Sect. [573] 


Bmean geometry 

Energy 1 

dipole 2 

quad . 2 

oct . 2 
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axi . 1 
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(%) 
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(%) 
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(%) 
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HD 35296: 2008 data with reduced number of profiles to match the 2007 dataset 
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1 this is a fraction (in %) of the total magnetic energy available. 

2 this is a fraction (in %) of the respective poloidal or toroidal field energy. 

Listed also is the fraction of the poloidal or toroidal magnetic energy in the dipolar (l = 1), quadrupolar (Ji = 2), octupolar (l = 3) and 
higher order (l > 4) components as well as the fraction of energy stored in the axisymmetric component (m = 0). 


HD 35296. The second star, HD 29615, was observed to be 
strongly poloidal with 75if per cent of the magnetic en¬ 
ergy being held in this configuration. Of this, 58 per cent of 
the poloidal component is dipolar (£ = 1). This is similar 
to other stars such as HD 76151 (SpType: G3V ) (|Petit et 
aT7||2008 1 , r Bootis (SpTy pe: F6IV) (|Fares et al.||2009 l and 
£ Bootis (SpType: G8V) ( |Morgenthaler et al.||2012| ) which 
are also strongly poloidal with complex fields with predomi¬ 
nantly dipolar fields. Table[5]gives a full listing of the various 
components for both HD 35296 and HD 29615. 


6.5 Differential rotation 

Differential rotation was measured for HD 35296 and 
HD 29615. The differential rotation of HD 35296, using 
Stokes V, was measured for the 2008 data using the % 2 
minimization technique, as described in Sect. |5.3] Using the 
magnetic signatures, HD 35296 has an equatorial rotational 
velocity, Q eq , of 1.804 ± 0.005 radd -1 with rotational shear, 
AQ, of Q.22 + _° 0 ° 0 t radd 1 . Differential rotation was not ob¬ 
served on the 2007 data due to the limited number of profiles 
obtained. This Afl value is consistent with similar observa¬ 


tions of other F-type stars by Reiners (2006) and Arnmler- 


von Eiff & Reiners (20121; however, Reiners (12006 1 found no 


evidence of differential rotation using the Fourier transform 
method of line profile analysis on HD 35296. Another star 
that exhibits this discrepancy is the pre-main-sequence bi¬ 
nary star HD 155555 (SpType = G5IV+K0IV) when [Duin 


stone et al. (2008) measured differential rotation on both 


components while 


Arnmler-von Eiff & Reiners (20121 could 


not detect differential rotation. Reiners|'p006f conducted a 


direct comparison between the Fourier transform method 
and Doppler imaging using HD 307938 (R58) in IC 2602. 
Marsden et al. (20051 used DI to measure a shear of Afl = 
0.025 ± 0.015 radd -1 . Using the Fourier transform method, 
the threshold for solid-body rotation, qg/qi, is 1.76 where 
qi and qg are the first two zeros of the line profile’s Fourier 


polar spot and small shear (such as R58), the spot has more 
influence on qg / qi than small deviations from solid-body ro¬ 
tation. This does not answer the question why HD 35296’s 
differential rotation was measured using ZDI but not Fourier 


transform method. Reiners (20061 measured qg/qi = 1.75 
(for HD 35296), only marginally less than solid-body ro¬ 
tation, yet this work measured a rotational shear, AD, of 
0.221o\q2 radd -1 . One would expect this to give a qg/qi 
~1.50 to 1.60. This disparity in measured differential rota¬ 
tion could be real, or just reflects the fact that this work is 


based on Stokes V whereas the value of Reiners (20061 was 


derived from Stokes / data which often gives lower values, as 
shown by the results for HD 29615. Alternatively, the Stokes 
/ LSD profiles of HD 35296 were not sufficiently deformed 
by spots to produce a detailed map whereas their presence 
may have affected the measurement in the Fourier domain. 

The differential rotation measurement for HD 29615, 
using the magnetic features indicate an equatorial rota¬ 
tional rate n eq = 2.74 Fq'qI radd -1 and shear Afi = 0.48ig;i2 
radd -1 . This rotational shear is relatively large although 
Reinhold, Reiners & Basri 1(2013} have measured larger ro¬ 
tational shear on active KEPLER stars with temperatures 
commensurate with HD 29615 (see Fig. 15 of that work). In 
contrast, the spot features provide a different value, namely 
n eq = 2.687 q;o2 radd -1 and Afl of 0.077o(o3 radd -1 (see Fig. 


131. Donati, Collier Cameron & Petit (2003b) noted that in 


some early K-dwarf stars, the level of differential rotation 
measured from brightness features is usually lower when 
compared with the levels measured when using magnetic 
features. Their interpretation of this variation is that the 
brightness features and the magnetic features are anchored 
at different depths within the convection zone of the star. 
This could also be the case here, though the difference in Afl 
that we measure for HD 29615 is more extreme than pre¬ 
viously observed for other stars. An alternative explanation 


transform. Reiners (2006 1 argued that for a star with a large 


may be offered by the work of Korhonen & Elstner (20111. 
Using dynamo calculations, they suggest that large starspots 
do not necessarily follow the actual differential rotation of 
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Figure 13. Differential rotation using Stokes I and V for 
HD 29615, including the associated variation ellipses (as ex¬ 
plained in Sect. |53) . The brightness (spot) features produce an 
equatorial rotational velocity, Cl eq of 2.68 ’ [['([(i rad <1 ! with rota¬ 
tional shear, A12 of 0.073,) 03 radd -1 . Using the magnetic features 
observed, HD 29615 has an equatorial rotational velocity, Q eq of 
2.74 radd -1 with rotational shear, AS1 of 0.48;(g^2 radd -1 . 


the star, but have more solid-body like behaviour, whereas 
the true surface differential rotation is only recovered if small 
magnetic features are added to the simulations. In addition, 
the paucity of low level spot features could affect the differ¬ 
ential rotation measurements. The differing measurements 
of differential rotation should be treated with caution until 
further Stokes I and V data are obtained for this star. 


6.6 Surface differential rotation, rotation rate and 
convection zone depth 


Barnes et al. (2005) showed a relationship between differen¬ 


tial rotation (Afi) and temperature, with increasing shear 
as a result of increasing temperature. Recent studies of G- 
dwarf stars such as HD 171488 (Marsden et al.|2006 Jeffers 


fc Donati|20 08 Jeff ers et al.|2011| ) and HD 141943 ( |Marsden 
et al. 2011b I show significant levels of differential rotation 


beyond those observed by Barnes et al. (20051. HD 29615, 
like HD 171488, has extreme levels of rotational shear with 
a = 0.48io!i2 radd -1 . HD 171488 has a 1.33 d rotation 
period with a usim = 38 kms -1 whereas HD 29615 has a 
rotational period of 2.34 d with ausini = 19.6 kms -1 . Given 
that both HD 171488 and HD 29615 have similar photo- 
spheric temperatures and age, it appears that rotation rate 
has limited effect on rotational shear, thereby supporting 


the observations by Barnes et al. (2005). 


Kiiker et al. ([20111, using theoretical models, demon¬ 


strated that the extreme surface shear of stars such as 
HD 171488 (and by implications HD 29615) can only be 
explained with a shallow convection zone. HD 171488 has a 


convection zone depth (CZD) of 0.206 R* (0.233 Rq) ( Jeffers 
et al.|2011) yet the CZD for HD 29615 was estimated to be 


0.252 ± 0.011 R*, as determined from the stellar evolution 
models of Siess et al. (2000 6 The estimate for HD 29615 
was based on the V-I colour index ( van Leeuwen|20'07 1 and 
the absolute magnitude, as determined from the maximum 
visual magnitude measured by the HIPPARCOS space mis¬ 
sion. O’Neal, Saar & Neff (19961 concluded that on some 


heavily spotted stars, the observed maximum V magnitude 
underestimates the brightness of the unspotted star by ~ 
0.3-0.4 mag. For example, HD 106506 ( |Waite et al. ||20 11 a 1 
had a visual magnitude of ~8.54 but the unspotted magni¬ 
tude was determined to be 8.38 (using the imaging code). As 
there was limited knowledge of the fractional spottedness of 
many of the stars that had their differential rotation already 
determined by other authors it was decided to use the maxi¬ 
mum visual magnitude listed in the HIPPARCOS database 
with the caveat that this is possibly leads to an overestimate 
of the true CZD. This measurement for HD 29615 places ad¬ 
ditional constraints on the conclusion of the theoretical work 
of Kfiker et al. (2011) that extreme levels of differential ro¬ 
tation require a shallow convection zone. 


7 CONCLUSIONS 

From the results presented in this paper, HD 35296 and 
HD 29615 are both Sun-like stars whose moderately rapid 
rotation has led to very complex surface magnetic fields that 
exhibit high levels of chromospheric activity and surface dif¬ 
ferential rotation. This variation in the chromospheric activ¬ 
ity, with rotational phase, was evidenced by modulation of 
the Ca 11 H & K, Ca 11 Infrared Triplet and Ha spectral 
lines (HD 35296) and Ha emission (HD 29615). High lev¬ 
els of differential rotation were measured. The differential 
rotation on HD 29615 showed a significant discrepancy in 
shear values between spot and magnetic features. This is an 
extreme example of a variation also observed for other lower- 
mass stars. There are indications that the magnetic field of 
HD 35296 appears to have undergone a minor reorganisa¬ 
tion of its field from 2007 to 2008 thereby demonstrating 
the evolving nature of the magnetic topology on this star. 
We conclude that the dynamo operating on both of these 
stars is similar to that of other more active, rapidly rotat¬ 
ing stars and is most likely a distributed dynamo operating 
throughout the convection zone. 
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